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The forbidden high ionisation line region of the type 2 
quasar Q1131+16: a clear view of the inner face of the 
torus? 
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ABSTRACT 

We present spectroscopic observations of the type 2 quasar SDSS J11311. 05+162739. 5 
(Q1131+16 hereafter; z=0.1732), which has the richest spectrum of forbidden high 
ionisation hues (FHIL, e.g. [Fe vii], [Fe x], [Fe Xl] and [Ne v]) yet reported for an 
AGN, as well as unusually strong [O iii]A4363 emission. The study of this object 
provides a rare opportunity to investigate the physical conditions and kinematics of the 
region(s) emitting the FHlLs. By comparison with photoionisation model results, we 
find that the FHIL region has high densities (10^'^ < uh < 10^ ° cm"'^) and ionisation 
parameters (-1.5 < log[U] < 0), yet its kinematics are similar to those of the low 
ionisation emission line region detected in the same object (FWHM ^ 360±30 km/s), 
with no evidence for a significant shift between the velocity centroid of the FHlLs and 
the rest frame of the host galaxy. The deduced physical conditions lie between those 
of the Broad-Line (nij>10^ cm"'^) and Narrow-Line Regions (n/f <10® cm"'^) of active 
galactic nuclei (AGN), and we demonstrate that the FHIL regions must be situated 
relatively close to the illuminating AGN (0.32 < rpHiL < 50pc). We suggest that the 
inner torus wall is the most likely location for the FHIL region, and that the unusual 
strength of the FHILs in this object is due to a specific viewing angle of the far wall 
of the torus, coupled with a lack of dust on larger scales that might otherwise obscure 
our view of the torus. 
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1 INTRODUCTION 

Most Seyfert galaxies show some spectral lines from forbid- 
den transitions of highly ionised ions in the ir spectra, e.g. 
[Fe v ii] , [Fe x], [Fe xi] and even [Fe xiv] jPenston et al.l 
11984 ). These emission lines are often blueshifted with re- 
spect to the rest frame of the AGN and have velocity 
widths which are between those of the Narro w-Line Region 
(NLR), and the Broad-Line Region (BLR. [ Penston et al.l 
1 1984 lAppenzeller fc WagneJ Il99ll . iMullanev et all l2009h . 
In some rare cases many forbidden high ionisation lines 
(FHILs) of relatively high equivalent w idth have been de - 
tected. Exampl es include III Zw 77 JOsterbrockl Il98lll . 
Tololo 0109-383 jFosburv fc Sansomiri983l ') and ESQ 138 Gl 
jAUoin et al.lll99a V 

The physical mechanisms and conditions that allow 
strong FHILs to be produced have been debated for some 
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time: whether there is a continuity from the photoionisa- 
tion proc esses which produce the lower ionisation line s in 
the NLR l|Korista fc Ferlandlll989l . iFerguson et al.lll99'7l ). or 
there is an entirely different mechanism for their forma- 
tion (e.g. coUisional excitation i n a high temperature gas; 
iNussbaumer fc Osterbrocklll970l '). 

Given that the transitions associated with many of the 
FHILs have high critical densities (nc>10^ cm~^), it has 
been suggested that they may originate in the innermost 
torus wall facing the illuminating source. Studies of the high 
critical density [Fe vil]A6086 emission line across Seyfert 
galaxies from types 1-2, including intermediate types, find 
that its strength increases relative to the low ionisation lines 
from type 2 to type 1. This can be interpreted in terms of 
the or ientation-based division between Seyfert types (e.g. 
lAnton ucci, 1993 ): as the Seyfert type gets closer to a Seyfert 
1, more of the emission fr om the inner torus becomes 
visible to the obs erver (see iMuravama fc Taniguchilll998l . 
iNagao et al.ll20()ll ). Therefore the high critical density lines 
that are preferentially emitted by the torus are likely to be 
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stronger in type 1 objects, as observed. However, to date the 
idea that the [Fe Vll] emission hnes are associated with the 
torus has not been thoroughly tested using detailed emission 
line ratios that measure physical conditions accurately. 

Although it seems plausible that at least some of the 
FHIL emission arises in the torus, the causes of the diversity 
in both the relative strength and kinematics of the FHILs re- 
main uncertain. One possibility is that the unusually strong 
FHILs observed in some objects are symptomatic of a recent 
energetic occurrence whi ch is able to illumina te the region 
that produces these lines (jKomossa et al.|[2009l ) ; suggestions 
for such events include the rapid accretion of material from 
the ISM, the tidal disruption of a stars by the supermas- 
sive black hole (SMB H), gamma-ray bursts and supernovae 
l|Komossa et al.ll2009l '). 

The unusual strength of the FHILs in objects such as 
HI Zw 77, Tololo 0109-383 and ESO 138 Gl provides us 
with a rare opportunity to study the nature of the region 
of the AGN which emits them, and thus helps us to better 
understand the structure of AGN in general. 

This paper reports an investigation of the object 
Q1131-I-16, which has the richest spectrum of FHILs yet re- 
ported for an AGN. Q11314- 16 was discovered in the 2MASS 
survey and was classified as a Seyfert 2 galaxy with a red- 
shift of 0.174. It belongs to sample of quasar-like objects 
identified on the basis o f their red near-IR colours (J-K>2.0) 
in the 2MASS survey (|Cutri et al.lll995l ). In this paper we 
present deep optical and infrared spectra of Q1131-I-16, and 
use these to deduce the physical conditions and kinematics 
of the FHIL emission lines, with the aim of investigating 
their origin. The cosmological parameters used throughout 
this paper are adopted f rom WMAP: Hp = 71 km s~^, 
= 0.27 and Qa = 0.73 (|Spergel et all 120031 '). resulting in a 
spatial scale of 3.46 kpc arcsec"^. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 ING ISIS observations 

Low resolution optical spectroscopic observations of 
Q1131-I-16 were taken on 9th February 2007 with the ISIS 
dual-arm spectrograph on the 4.2-m William Herschel Tele- 
scope (WHT) on La Palma as part of a spectroscopic sur- 
vey of a representative sample red, quasar-like objects. The 
full sample for this survey comprised a complete RA-limited 
sub-sample of 2 4 objects with (J-K>2 ) and z < 0.2 selected 
from the list of lHutchings et al]|2003l . which is itself repre- 
sentative of the population of red, 2MASS-selected quasars. 
The unusually strong FHIL in Q1131-I-16 were discovered 
serendipitously in the course of this survey. 

In the red, the R158R grating was used with the RED- 
PLUS CCD, and in the blue, the R300B grating was used 
with the EEV12 CCD. A dichroic at 5300A was employed 
to obtain spectra with useful wavelength ranges ~3250- 
5250A in the blue, and ~5200-9500A in the red. To reduce 
the effects of differential refraction, all exposures were taken 
when Q1131-I-16 was at low airmass (sec z < 1.1) and with 
the slit aligned close to the parallactic angle. The seeing for 



the night of the observations varied from 0.8 to 1.3 arcsec- 
onds (FWHMfl 

Sets of three 600s exposures were taken on both arms 
simultaneously, giving a total exposure time of 1800s. The 
data were taken with a 1.5 arcsecond slit along a position 
angle PA315. To eliminate contamination from second order 
emission, a GG495 blocking filter was introduced into the 
ISIS red arm. 

The data were reduced in the standard way (bias sub- 
traction, flat fielding, cosmic ray removal, wavelength cali- 
bration, flux calibration) using packages in iRAlQ. The two- 
dimensional spectra were also corrected for spatial distor- 
tions of the CCD. To reduce wavelength calibration errors 
due to flexure of the telescope and instrument, arc spec- 
tra were taken at the position of the object on the sky; the 
estimated wavelength calibration accuracy is O.lA in both 
the blue and red respectively (however, this may be an un- 
derestimate at the extreme edges of the spectra). The atmo- 
spheric absorption features were removed by dividing by the 
spectrum of a telluric standard (BD+17 2352), taken close 
in time and airmass to the observations of Q1131-I-16. The 
spectral resolution, calculated using the widths of the night 
sky emission lines, was 6.11 ± 0.74A in the blue and 10.96 
± 0.83A in the red (measured in the observed frame of the 
spectrum). The spatial pixel scales of the 2D spectra are 0.4 
arcseconds in the blue and 0.44 arcseconds in the red, and 
the relative flux calibration uncertainty - based on 11 ob- 
servations of 8 flux standard stars taken throughout the run 
- is estimated to be ±5%. 

The spectra were extracted and analysed using the 
STARLINK packages FIGARO and dipso. 



2.2 Gemini GMOS observations 

In order to test the possibility that the emission lines of 
Q1131-I-16 are variable, a follow-up spectrum was taken us- 
ing GMOS on the Gemini south telescope at Cerro Pachon. 

Spectroscopic and imaging observations using GMOS 
were taken on the 21st February 2010 during Gemini queue 
observation time as part of the program GS-2009B-Q-87. 
The spectral observations comprised three 600 second ex- 
posures using a central wavelength 5000A at an airmass of 
1.46, and a 1.5 arcseconds slit aligned along position angle 
PA163° — close to the parallactic angle for the centre of 
the observations. This resulted in a spectrum with a use- 
ful wavelength range ~3700-6450A. The observations of the 
flux standard star, Feige56, had the same set-up but with 
an exposure time of 10 seconds, taken at an airmass of 1.61 
and with a 5 arcsecond slit. Deep GMOS images were also 
taken using an r' filter at an airmass of 1.56. Each exposure 
of the imaging observations had an exposure time of 250 sec- 
onds, and the data were taken in a 4 point dither pattern, 
resulting in a total exposure time of 1000s. The seeing was 



^ Unfortunately, due to a lack of suitable imaging observations, 
an accurate seeing estimate does not exist for the exact time of 
the WHT observations of Q1131+16. 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
the Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation ( |http:/ /iraf.noao.edu/p . 
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estimated to be 0.75±0.02 arcseconds (FWHM), based on 
measurements of stars in the GMOS image. 

The spectra were extracted and analysed using the 
STARLINK packages FIGARO and dipso, and the IRAF pack- 
ages GEMINI and NOAO. Analysis of the night sky lines im- 
plies a wavelength calibration uncertainty of 0.50±0.06A for 
the spectrum. The spectral resolution, calculated using the 
widths of the night sky emission lines, was 7.12 ± O.llA m 
the observed frame of the observations, and the spatial pixel 
scale is 0.147 arcseconds per pixel. 

Comparisons between the WHT and Gemini spectra 
demonstrate that the fluxes of the stronger emission lines 
(e.g. [O iii]AAA4363,4959&5007, [Nev]A3426, [Fe vii]A 3759) 
and continuum agree within 20% between the two sets of 
observations. This is remarkable considering that the two 
sets of observations were taken with different instruments, 
at different slit position angles and under different seeing 
conditions. 

The GMOS image was reduced using the standard Gem- 
ini pipeline reduction software. The photometric calibration 
used the zero point magnitude for the r' filter advertised on 
the GMOS south website, taking full account of the airmass 
and exposure time of the observations. 

2.3 ING LIRIS observations 

On the night starting the 2nd March 2010 we obtained a 
1.388 - 2. 419/jm spectrum of Q1 131-t-16 using the HK grism 
in LIRIS l|Manchado et al.ll 19981 ') on the 4.2-m William Her- 
schel Telescope on La Palma. The spectrum has a dispersion 
of 9.7A per pixel, and a 1 arcsecond slit was used, giving a 
spectral resolution of ~40A. The seeing was approximately 
2 arcseconds, and the night suffered from heavy cirrus. A to- 
tal exposure time of 88 minutes was obtained over airmasses 
ranging from 1.0 to 2.2. The total exposure was divided into 
240 second sub-exposures and the telescope was nodded in 
the standard manner to aid sky subtraction. A random jit- 
ter of 10 arcseconds was added to the default nod positions 
to limit the impact of bad pixels. Observations of the A5V 
star BD-l-16 2325 were also taken to correct for the effects of 
telluric absorption and to provide a relative flux calibration. 
Both stars were observed with a slit position angle of zero 
degrees. Fo r a more detailed descriptio n of the reduction 
process, see IKamos Almeida etHI l|2009l) . 



3 RESULTS 



3.1 Spatial distribution of the emission lines 

While our 2D spectra show that all the emission lines are 
strongly concentrated on the nucleus of the host galaxy, we 
detect spatially extended line emission in [O ll]A3727 (see 
Figure 1), H/?, [O iii]AA5007,4959 and the Ha+[N ii] blend. 
The emission lines have a maximum extent of 9 arcseconds 
in the north-south direction (29kpc) in the Gemini spec- 
trum (PA163), and 15 arcseconds in the north west-south 
east direction (47kpc) in the WHT spectra (PA315), corre- 
sponding to the spatial extent of the galaxy disk visible in 
our Gemini r' image (see ^3.5). 



In order to investigate the nuclear spectrum, we ex- 
tracted apertures of size of 1.5x1.5 arcseconds from the 2D 
WHT and Gemini spectra. The apertures were centred on 
the nucleus. 

To determine the spatial distributions of both the emis- 
sion lines and continuum along the slit, spatial slices were 
extracted from the 2D WHT spectrum over the wavelength 
ranges given in Table 1. The continuum slices were extracted 
with similar wavelength ranges to their nearby emission lines 
so that the continuum could be accurately subtracted from 
the slices containing the emission lines. The dipso STARLINK 
package was then used to measure the centroid and spatial 
FWHM of the flux distribution of each slice, by fitting a 
single Gaussian profile to the central cores of the emission. 

Overall, the centroids measured for each emission line 
and continuum slice in Table 1 are all consistent within their 
uncertainties, implying no spatial offset between any of the 
emission lines individually, nor between the emission lines 
and the continuum. 

The spatial FWHM of all the emission lines — with the 
exception of [O ll]A3727 — are consistent with the seeing 
on the night of the observations (see ^2.1), suggesting that 
their spatial distributions are unresolved in the observations. 
Therefore the spatial distributions of the high ionisation 
lines provide an indication of the true seeing of the WHT ob- 
servations (1.1 < FWHM < 1.35 arcseconds). As expected, 
there is evidence that the seeing degrades towards shorter 
wavelengths, since the two shortest wavelength FHIL have 
significantly larger spatial FWHM than their longer wave- 
length counterparts. Moreover, fact that the spatial FWHM 
measured for all continuum slices significantly exceeds the 
estimated seeing, demonstrates that the continuum emission 
is spatially resolved in the nuclear regions. 

The spatial FWHM of the [O ll]A3727 emission line is 
significantly broader than the other FHILs at similar wave- 
lengths in Table 1, and this line is clearly resolved in the 2D 
spectrum (see Figure 1). The combination of the large spa- 
tial extent and spatially resolved nuclear emission of [O ll] 
makes this emission line ideal for the study of the emission 
line kinematics across the host galaxy Q1131-I-16 presented 
in ^3.6. 

3.2 Line identifications 

The optical spectra of Q1131-I-16 are shown in Figures 2- 
5. What makes the spectrum of Q1131-I-16 special is not 
only the large number of emission lines, but also their 
rich variety. In particular, the detected FHILs include 
[Fe v]AAA3839,3891&4181, [Fe vi]AAA5146,5176&5335, [Fe 
vii]AAAAAA3759,4893,5159,5276,5720&6086, [Fe x]A6375, 
[Fe xi]A7891 and [Ne v] AA3346&34243. As well as the FHILs, 
lower ionisation species which are not typically detected 
in quasar spectra are found in the spectrum, including [Fe 
lv]AAAA2829,2836,4903&5236, along with Bowen resonance 
fluorescence lines of O ill, such as O ill A3133. Also de- 
tected are more typical AGN emission lines such as Ha, H/3, 
He II A4686, [O ii]AA3726,3729, [O iii] AAA4363,4959&5007. 

^ In this study we define a FHIL as an emission species with an 
ionisation potential greater than or equal to 54.4 eV (that of He 
II). 
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Table 1. Centroid positions and spatial extents of selected emission lines and the continuum at different wavelength intervals, as measured 
from the 2D WHT spectrum. The centroids (measured in pixels) and the FWHM (measured in arc seconds) were determined using single 
Gaussian fits to the bright cores of the s])atial <listril)utioiis of the omission linos along tlio slit. 





Range A 


Centroid 


± 


FWHM (arc seconds) 


± 


Blue Arm 


[Ne V] 3425 


35 


98.34 


0.01 


1.35 


0.01 


cont. 3392 


35 


98.23 


0.04 


1.71 


0.04 


[O 11] 3727 


45 


98.39 


0.08 


1.56 


0.06 


cont. 3690 


45 


98.21 


0.07 


1.78 


0.06 


[Fc VII] 3759 


45 


98.32 


0.07 


1.28 


0.05 


[Ne 111] 3868 


45 


98.32 


0.01 


1.24 


0.01 


cont. 3923 


45 


98.33 


0.06 


1.81 


0.07 


H5 


50 


98.49 


0.05 


1.04 


0.05 


cont. 4147 


50 


98.25 


0.06 


1.74 


0.06 


H7 


35 


98.33 


0.04 


1.09 


0.04 


[O 111] 4363 


35 


98.35 


0.05 


1.14 


0.02 


cont. 4466 


35 


98.38 


0.05 


1.73 


0.07 


Red Arm 


He II 4686 


50 


89.82 


0.12 


1.25 


0.15 


cont. 5545 


50 


89.84 


0.04 


1.74 


0.05 




50 


89.74 


0.04 


1.09 


0.04 


[O 111] 4959&5007 


130 


89.75 


0.01 


1.18 


0.01 


cont. 5093 


130 


89.85 


0.04 


1.74 


0.05 


[Fe VII] 6086 


65 


89.83 


0.02 


1.19 


0.02 


cont. 6139 


65 


89.80 


0.04 


1.87 


0.05 


Ha 


145 


89.81 


0.01 


1.09 


0.02 


cont. 6435 


145 


89.81 


0.04 


1.71 


0.05 


[Fe XI] 7892 


40 


89.76 


0.10 


1.18 


0.11 


cont. 7927 


40 


89.84 


0.06 


1.80 


0.07 




NW 



Figure 1. A section of the 2D WHT spectrum along PA315 showing [O n]A3727 and [Fe Vll]A3759. Note that the [O 11] is extended 6.9 
arcseconds above and 8.1 arseconds below, the continuum of Q1131+16 and is well resolved. The spatial direction of the slit is indicated 
along side the plot. In addition, a nuclear velocity gradient in the [O 11] can be seen from the fact that its profile is slanted when compared 
to that of [Fe vii] 
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However, the FHILs are unusually strong with respect to 
these latter lines, with [Ne v] A3426 comfortably exceeding 
the strength of the [O n]A3727 doublet, and [Fe vn]A6086 
of comparable strength to H/3. Also notable is the unusual 
strength of [O lll]A4363 compared to H7 and [O lll]A5007 
([O iii]5007/4363 = 5.86±0.21 in the WHT spectrum), as 
well as the relatively modest ratio of [O lll]A5007 to H/3 ([O 
ni]/H/3=5.27±0.16), given the strength of the other high 
ionisation lines. 

The strongest emission lines in this object have rela- 
tively large equivalent widths (EWs), including a large num- 
ber of FHILs with EWs comparable to other typical lower 
ionisation emission lines. Indeed, the EWs of the FHILs in 
Q1131+16 are larger than those of any other AGN with pub- 
lished spectra, including other objects with unusually strong 
FHI Ls such as III Zw 77 (IOsterbrocklll98ll 'l and Tololo 0109- 
383 (iF osburv fc Sanso ml ll983l ). Moreover the [O ill] emission 
line luminosity of Q1131-I-16 (4.6±0.1xl0* Lq) would lead 
to its class ification as a quasar 2 object according to the 
criterion of lZakamska et al] (|2003l ). Overall, the rich variety 
of emission lines, and their relatively large EWs, allow for 
a thorough investigation of the physical conditions of the 
FHIL emission region. 

A full list of line identifications made from these spec- 
tra is presented in Table fQ AH the line identifications have 
been determined by fitting single Gaussians to the emission 
features in both the WHT and Gemini spectra. The rest- 
frame wavelength ranges that were fitted for Q1131-I-16 are 
2700-4600 A and 4400-8000 A in the WHT spectra, and 
3100-5450 A in the Gemini spectrum. Table 5 gives the line 
flux ratios relative to H/3 for both spectra. The fluxes of the 
emission lines have not been corrected for intrinsic redden- 
ing for reasons that will become clear in /3.8. In addition, 
it appears that there is no need of a Galactic extinction cor- 
rection, since the IRSA extinction tool in the NED giv es a 
reddening of only E(B-V)= 0.0306 jSchlegel et al.lll998l ). 

The majority of line identifications have been confirmed 
in other astrop hysical objects such as AGN and plane- 
tary nebula (e.g. Osterbrock'l98ll. lFosburv fc Sanso mi l 19831 . 
lAlloin et al.] ri992. Kalcr 1976). However, there are several 
emission lines on Table 5 which have been identified using 
the National Institute of Standards and Technology (NIST) 
spectral line database (see ref. 4 in Table 5). An emission 
line was only regarded as a secure ID if its line centre was 
within 1.5 sigma of the wavelength predicted for that par- 
ticular emission line based on the mean redshift, and if its 
S/N ratio exceeded 3.0. As well as the 66 identified emission 
lines, there are 37 emission lines which remain unidentified 
(see Table 6), in the sense that there are no IDs for them in 
the NIST spectral line database that give redshifts within 
1.5 sigma of the mean redshift, yet their S/N exceeds 3. Note 
that, due to its higher S/N, the Gemini spectrum reveals a 
large number of faint lines that were not detected in the 
original WHT spectrum (see Tables 5 and 6). 

Interestingly, although there are many iron FHILs in 
the spectrum of Q1131+16, there is no clear evidence for 
the [Fe xiv]A5303 emission line. The identification of this 
emission feature at ~5300 A has been controversial in the 



past: it has b een debated whether i t is [Fe xiv1A530 3 or [Ca 
v]A5309 (e.g. lOke fc SargentllTgei . IWeedmanlll97lh . When 
this feature is fitted in our WHT spectrum, if only the [Ca 
v]A5309 identification is considered, the individual redshift 
of the emission line (0.17319±0.00011 for the WHT spec- 
trum) agrees within the uncertainties with the average red- 
shift of aU the fines (0.17325±0.00001, see /3.6). However, 
if the feature is identified with [Fe xiv]A5303, the individ- 
ual redshift becomes 0.17455±0. 00011, which is significantly 
higher (>10(t) than the mean redshift. 

A further interesting feature of the spectrum is that 
the Ha-|-[N 11] emission blend shows tentative evidence 
for a broad base of rest- frame width 11,500±2200 km s~^ 
(FWHM, see Figures 4 fc 6). This is consistent w ith the 
presence of a scattered (|Antonucci fc Milled [l985l) . or di- 
rectly observed, broad-line region (BLR) component. Spec- 
tropolarimetry observations will be required to confirm the 
scattered BLR possibility. 



3.3 Infrared spectrum 

In order to further understand the nature of this object we 
made spectroscopic observations of Q113H-16 at near-IR 
wavelengths. The observed wavelength range (1.4 to 2.4/im) 
was selected in order to simultaneously detect both Paa and 
Pa^. 

The K-band NIR spectrum of Q1131+16 taken using 
LIRIS is presented in Figure 7. Due to its relatively low 
S/N, this spectrum does not show the abundance of emis- 
sion lines seen at optical wavelengths. Although weak nar- 
row Paof and Pa4f| emission lines are detected, there is no 
sign of any broad Paa and Pa/? components, implying the 
BLR is enshrouded by dust. A single Gaussian fit to the 
PaQ line shows that it is unresolved within the uncertain- 
ties for the resolution of our observations, and its redshift 
(z=0.17314±0. 00011) is consistent with those of the optical 
emission lines. 

Any quasar nucleus component present in this sys- 
tem must be highly extinguished; the relatively red near-IR 
colours of this source measured by 2MASS do not appear to 
be due to a moderately extinguished quasar component that 
becomes visible at the longer near-IR wavelengths. There- 
fore this object cannot truly be described as a 'red quasar'. 



3.4 Spectral fitting model 

The emission lines in the spectrum of Q1131-I-16 were ini- 
tially fitted with single Gaussian profiles, but such fits did 
not provide an entirely adequate fit to the wings of the 
stronger lines (see Figure 8). To overcome this, a double 
Gaussian model was fitted to the spectral lines. This model 
is based on the fit of the [O iii]AA5007,4959 lines, because 
they are particularly strong emission features. To produce 
the model, both broad and narrow Gaussians were fitted to 
the Gemini [O ill] line profiles, where the centres, widths 
and intensities of each Gaussian were free parameters. The 
widths (FWHM) of the narrow components for the [O ill] 



* Unfortunately tables 5 & 6 are not available here, however they 
are available in the source tar. 



^ The Pa/3 emission line is on the blue edge of the H-band NIR 
spectrum, because of this any measured emission line data is ques- 
tionable. 
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Figure 2. WHT nuclear spectrum of Q1131+16 taken on the ISIS blue arm. Double lines which are labeled with an indicates that 
both components are resolvable, those lableled with a '/' indicates that both components are not resolvable. The flux scale is measured 
in units of lO"'^® ergs cm~^. Note the strength of high ionisation lines such as [Ne v] and [Fe vii] relative to [O ii], and of 

[O iii]A4363 relative to H5. For reference we also show a scaled version of the night sky spectrum extracted from the 2D frames at the 
bottom of the plot. 




Figure 3. WHT nuclear spectrum of Q1131+16 taken on the ISIS red arm. Double lines which are labeled with an '&' indicates that 
both components are resolvable, those lableled with a '/' indicates that both components are not resolvable. The flux scale is measured 



Q1131+16: a clear view of the inner face of the torus? 7 

Q1131 + 16 - Aperture 1 




4800 5000 5200 5400 5600 5800 6000 6200 6400 6600 6800 7000 7200 7400 7600 7800 



Rest ¥avelength(A) 



Figure 4. An expanded plot of the Q1131+16 ISIS red arm spectrum. This has been presented to highlight the weaJcer features in the 
red spectrum. The flux scale is measured in units of 10~^® ergs cm~^. As well as the multitude of high ionisation lines, note 

the broad base to the Ha+[N ii] blend. For reference we also show a scaled version of the night sky spectrum extracted from the 2D 
frames at the bottom of the plot. 
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Figure 5. Nuclear spectrum of Q1131+16 taken using GMOS on the Gemini South telescope. The flux scale has been narrowed here to 
highlight how rich the spectrum is in emission lines. The flux scale is measured in units of lO"'^^ ergs A~^ cm~-^. The gaps in the 
spectrum are due to the gaps between the CCD chips in the GMOS instrument. 



M. Rose, C. N. Tadhunter, J. Holt, C. Ramos Almeida and S. P. Littlefair 



J 




Figure 6. The broad component centered on the Ha+[N ii] emission blend. The telluric feature is fitted as well as the various emission 
components from Q1131+16. 
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Figure 7. The K-band LIRIS spectrum of Q1131+16. The rest wavelength range shown here is 17000-20500 A. The only significant 
feature is the Paa emission line. 
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Figure 8. Fits to the [O III] AA5007, 4959 emission lines. Top. Single Gaussian fit to both [O in]A5007 and A4959. The overall fit does not 
fit the wings of the emission lines well. Bottom. The double Gaussian fit to the [O in]A5007 and A4959 emission lines from the Gemini 
data. Here the model has had more suecess in fitting the wings of the emission lines when eoniparcd to that above. The individual 
components of the fit are included on the plot; the components of A5007 are drawn with a dotted line and the components of A4959 are 
fitted with a dashed and dotted line. The overall line profile is fit with a dashed line. The model throughout this investigation is based 
on the [O in]A5007 emission line. The flux is measured in 10~^^ ergs s~^ A"'^ cm~-^ and the wavelength is in units of A. 
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emission lines were found to be consistent with the instru- 
mental width of the Gemini spectra (7.1 A) at 6.9±0.2 A. 
We therefore use the instrumental width to fit the narrow 
components in the model fit to other lines, varying this to 
take into account the resolution of the different spectra. 
The measured FWHM of the broad component (14.1±0.6 
A(FWHM), corresponding to a rest-frame velocity width of 
720±30 km s~^) was used to obtain an intrinsic velocity 
width for the broader component by correcting its FWHM 
in quadrature using the instrumental width. In addition, the 
broad component is redshifted by 1.8±0.1 A from the nar- 
row component, corresponding to a velocity shift of 92±6 
km in the galaxy rest frame. The double Gaussian [O 
in] model was then fitted to all other emission lines, in order 
to derive the line fluxes listed in Tables 4 and 5. 

Figure 8 shows the spectral fit to the [O in]A5007 emis- 
sion line used to determine the spectral fitting model. The 
overall fit in the bottom panel of Figure 8 is significantly 
better than that shown in the top panel, where the emission 
lines are fitted with a single component free fit. The param- 
eters from the double Gaussian [O ill] fits were used to fit 
the emission lines throught the spectrum. In general these 
fits were extremely successful. There is, however, a minor- 
ity of lines — indicated by notes 1 and 2 in Tables 4 and 5 
— which are not fitted well by both components of the [O 
III] model. These lines are generally weak, and many are in 
blends with other emission lines. 

The [N ii]AA6548,6584 doublet is blended with the Hq 
emission line and therefore was modelled using constraints 
provided by atomic physics (i.e. same FWHM, 1:3 intensity 
ratio, and line centre of A6548 fixed to the centre of A6583). 
This is true of other doublet blends, for example the [O 
i]AA6300,6364 (blended with [S iii]A63i2 and [Fe x]A6375). 

3.5 Host galaxy morphology 

Figure 9 presents our Gemini r'-band image of Q1131-I-16. 
The galaxy morphology appears consistent with a late type 
spiral classification (Sb or Sc), with a moderate degree of 
asymmetry in the sense that the spiral arms are more ex- 
tended to the west than the east of the nucleus. However, 
we do not find evidence for the "large scale and clear tidal 
disturbance" reported bv lHutchings et "all (|2003h . Clearly, a 
more detailed analysis of the morphology and kinematics of 
the host galaxy is required in order to determine the true 
interaction status of this system. 

Based on aperture photometry, the total magnitude of 
the host galaxy was estimated using a circular aperture of di- 
ameter 23 arcseconds (57 kpc), giving an r'-band magnitude 
of 17.11±0.26 — consistent with the Petrosian r-band mag- 
nitude for this object listed on the SDSS web site. Note that 
this estimate represents an upper limit on the brightness of 
the stellar component of the host galaxy since no correction 
has been made for AGN components such as emission lines 
and scattered light. 

3.6 Emission line kinematics and redshift 

The median redshift determined using all of the securely 
identified emission lines is z=0.I7325±0. 00001. From only 
the securely identified FHILs (IP ^ 54.4 eV) the redshift is 



z=0.17323±0. 00002, while for the low ionisation lines we find 
z=0.17326±0. 00001. These values are consistent within their 
uncertainties, as emphasised by Figure 10, where estimated 
redshifts of all the lines are plotted. 

In order to check whether the mean emission line red- 
shift represents that of the host, the [O ll]A3727 emis- 
sion line blend was fitted as a function of position across 
the galaxy for both Gemini and WHT spectra. Figure 11 
presents the results. The velocities were determined relative 
to the median redshift of all the emission lines of the system 
(z=0.17325±0. 00001), assuming a rest-frame wavelength of 
3727.4 A for the blend — representing the flux weighted 
wavelength of the [O ll] AA3726,3729 blend for a density of 
10"^ cm~^. For both position angles, the central region (±1.2 
arcseconds, ±4kpc) shows a velocity field that resembles a 
rotation curve; the velocity gradient of the [O ll]A3727 emis- 
sion in the nucleus can also clearly be seen in Figure 1. This 
component of emission may correspond to an inner disk of 
the galaxy, implying that the [O ll]A3727 emission compo- 
nent in Q1131-I-16 is resolved in the nucleus. However, the 
two datasets differ significantly at larger radii; such differ- 
ences are likely due to the different slit PAs, given that each 
slit intersects different regions of Q1131-I-16 see Figure 11): 
whereas the slit for the Gemini data is within ~ 20 degress 
of the minor axis of the host galaxy, that for the WHT is 
offset for the minor axis by ~ 50 degress. 

The rotation curve dereived from the PA315 WHT data 
appears flat at large radii. However, the velocities of the 
flat parts of the velocity curves are asymmetric about the 
host redshift measured from the nuclear spectrum: the rota- 
tion velocity ranges from -|-90±5 km/s for the region mov- 
ing away from the observer to -176±5 km/s for the region 
moving towards the observer in the WHT spectrum. This 
asymmetry can be explained by either the uncertainty in the 
electron density assumed in determining the rest wavelength 
of the [O ll] blend (rie is likely to differ between the nucleus 
and the disk of Q1131+16), or the uncertainty in the wave- 
length calibration of the spectrum. Both of these aspects 
can systematically shift the curve in Figure 11, leading to 
an asymmetric appearence. However, despite the apparent 
asymmetry, the bisector of the two flat parts of the rotation 
curve is within 50 km s~^ of the rest frame of the host galaxy 
estimated on the basis of the median redshift derived from 
the nuclear spectrum. Thus our data provide no evidence for 
high velocity outflows in the nuclear emission line regions. 

On the other hand, the rotation curve derived from the 
PA163 Gemini data appears more complex: outside the in- 
ner rapidly rotating region, the velocities drop signiflcantly 
at intermediate radii (1.2<r<2.5 arcseconds) on either side 
of the nucleus, before rising again at larger radii. It is pos- 
sible that the inner rapid rotation along PA163 samples the 
seeing disk of the same inner disk component detected in 
the WHT data along PA315. In this case, the initial drop in 
the velocities (1.2<r<2.5 arcseconds) may be due to the slit 
sampling the outer disk close to the minor axis where the 
(projected) rotation velocity amplitude is likely to be small. 
However, given the complexity of the rotation curve along 
PA163, non-circular gas motions cannot be ruled out. 

Our double Gaussian [O lll] model fits the majority 
of the lines well in the nuclear spectrum, suggesting that 
the velocity widths of the FHIL and low ionisation emission 
lines are similar. The median rest-frame velocity width of 
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the emission lines, based on single Gaussian fits, is found to 
be 361±29 km s"^ (FWHM). 

Figures 12 and 13 also show that there is no correlation 
between the widths of the emission lines as measured by sin- 
gle Gaussian fits, and their critical densities and ionisation 
potentials. Again this suggests the kinematics of the FHILs 
and the low ionisation lines are similar. 



3.7 FHILs in 'typical' Seyfert galaxies 



ectra of all types of Seyfert 
These emission lines have 



FHILs are detected in the sp e 

galaxies (jPenston et al.l 1 1984 ). 

been studied rigorously in the past (e.g . Periston et al.l 
1 1984 iNagao et al]|20od iMullanev fc Wardll2008l ). and sev- 



eral common properties have been identified. It is therefore 
important to highlight the differences between the spectrum 
of Q1131+16 and the spectra of 'typical' Seyfert galaxies. 

Possibly the most notable property of the FHILs in typ- 
ical AGN is that the velocity widths (FWHM as measured 
by single Gaussian fits) tend to be intermediate between the 
NL R and BLR (in the range of 500 < FWHM < 1000 km 
s~^. lAppenzeller k, 6streicherlll99ll ). However the FHILs in 
Q1131-I-16 and other similar objects (e.g. HI Zw 77 and 
Tololo 0109-383 etc.) do not share this property. In the case 
of Q1131-I-16, Figures 12 and 13 show no significant varia- 
tion in the velocity width (FWHM) with the critical density, 
and ionisation potential. 

For typical AGN it is also found that FHILs are 
blueshifted with respect to the host galaxy rest-frames. 
In fact, the blueshift correlat es with the observed velocity 
FWHM (jMullanev et al.ll2009l ). There is no evidence for this 
in the spectrum of Q1131-I-16. The individual redshifts from 
single Gaussian fits of each emission line are plotted in Fig- 
ure 10. From this it is clear that there is no tendency for the 
FHILs to be blueshifted with respect to the average redshift 
of the low ionisation lines; the points representing the FHIL 
redshifts are scattered about the average redshift of the low 
ionisation species. 

Finally, the FHILs of Q113H-16 and similar objects 
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Figure 10. Individual redshifts of the emission lines of Q1131-I-16 
as measured for the WHT data. The triangles indicate the low 
ionisation emission lines and neutral species emission, the black 
circles indicate the FHILs. All points lie within 1.5 a of the 
mean redshift, however the longer wavelength emission lines ([O 
ii]AA7320&7330, [Ar iv]A7868 and [Fe xi]A 7892) do not. This 
is likely due to larger wavelength calibration errors at the long 
wavelength end of the spectrum. 



have larger equivalent widths when compared to those of 
other 'typical' Seyfert galaxies. 



3.8 Balmer decrements 

In order to derive accurate physical conditions, it is crucial to 
correct for any intrinsic dust extinction. The most common 
way to achieve this for optical spectra is to use Balmer line 
recombination ratios. 

Table 2 presents flux ratios for the Balmer recombina- 
tion lines measured from the Gemini and WHT spectra. The 
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Figure 11. The rotation curves for Q1131+16. The crosses 
(WHT, PA315) and filled squares (Gemini, PA163) are the mea- 
sured velocity shifts determined at different spatial locations 
along the slit of the long slit spectrum. The directions of the 
slits are indicated in the upper left-hand corner. Velocity shifts 
are determined by the difference between the measured redshift 
of [O ll] on each pixel and the average redshift of the emission 
lines. The horizontal line indicates the bisector of the flat parts 
of the rotation curves at large radii along PA315. 
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Figure 12. Rest-frame line width (FWHM) versus critical den- 
sity for the emission lines of Q1131-I-16. The line widths have 
been determined from single Gaussian free fits to the emission 
lines, and have been corrected for the instrumental profile. 
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Figure 13. Rest-frame line width (FWHM) versus ionisation po- 
tential for the emission lines of Q1131-I-16. The line widths have 
been determined from single Gaussian free fits to the emission 
lines, and have been corrected for the instrumental width. 



Table 2. Line flux ratios for Balmer recombination emission lines 
for both the WHT and Gemini observations. Ratios are relative to 
H/3. The last three columns on the table are the expected ratios 
for case B recombination at different tem peratures at a typical 
NLR density (10^ cm~ 3. IOsterbroc3ll978h . 



Species 


WHT 


Gemini 


lOOOOK 


20000K 


Ho 


4.99±0.17 




2.85 


2.74 


H7 


0.44±0.02 


0.45±0.03 


0.47 


0.48 


H<S 


0.24±0.02 


0.24±0.01 


0.26 


0.26 



H7/H/3 and H5/H/3 ratios in both spectra agree with their 
theoretical Case B values within the uncertainties, assum- 
ing conditions typical of the NLR. This suggests that there 
is no need for an extinction correction for the spectrum of 
Q1131-I-16. Note that, because of the large equivalent widths 
of the Balmer lines, uncertainties due to underlying Balmer 
absorption lines from young stellar populations in the host 
galaxy are not a serious issue for Q113H~16. 

Another interesting feature of the spectrum of 
Q1131+16 is that its Ha flux is high when compared to other 
Balmer lines in the series. The Ha/H/3 ratio is not consistent 
with Case B recombination (Ha/H/3=4.99±0.17), see / 4.2 
for a detailed discussion. Given the pote ntial for collisional 
excitation of the Ha line (see (j> 4.2) (jCaskell fc Ferlandl 
.1984). we ignore the Ha/H/3 ratio when assessing the level 
of extinction in Q1131-I-16. Therefore, since the high order 
Balmer line ratios for this source provide no evidence for 
significant reddening, we make no reddening correction of 
the line ratios. 
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Table 3. Diagnostic line intensity ratios measured for the nuclear 
regions Q1131+16, as derived from the WHT data. Values in 
brackets are the ratios given by the Gemini data. Potentially, the 
[Fe vil]A6086 line may suffer from some contamination by [Ca 
v]A6087. However, the fact that the [Fe vil] (5720/6086) ratio 
is within 1 c of the value predicted by atomic p hysics ([Fe VII 
(5720/6086) = 0.617: iNussbaumer fc Storevl[T983) indicates that 
any such contamination is modest. 
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Species 


Ratio 


± 


[O ii](3726+3729)/(7317+7330) 


5.02 


1.19 


[O III] (5007/4363) 


5.86(6.84) 


0.21(0.32) 


[O iii](5007/H/3) 


5.27(5.74) 


0.16(0.16) 


[S ii](6717+6731)/(4069+4076) 


7.51 


1.46 


[Fe vi](5176)/[Fo v](3891) 


0.267(0.176) 


0.076(0.008) 


[Fe vii] (6086/3759) 


1.18 


0.076 


[Fe vii] (5159/6086) 


0.208 


0.067 


[Fe vii] (5720/6086) 


0.667 


0.067 


[Fe x](6375)/[Fe vii](6086) 


0.473 


0.074 


[Fe xi](7892)/[Fe vii](6086) 


0.317 


0.046 


[Ne v](3426)/[Fe vii](6086) 


2.04 


0.207 
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Figure 14. Diagnostic plot showing the 

measured [Sll](6717+6731)/(4069+4076) and 

[Oll](3726+3729)/(7317+7330) transauroral ratios (filled point), 
compared with photoionisation model results for different densi- 
ties and ionisation parameters. Equal densities are connected by 
the solid lines (as labelled), and the ionisation parameters are 
connected by the dotted lines (also labelled). 



4 DISCUSSION 

4.1 Physical conditions implied by the FHILs 

To investigate the physical conditions implied by both the 
low ionisation and high ionis ation species, the ph otoionisa- 
tion model code CLOUD-^ (iFerland et al.lll998l ) was used 
to create single slab photoionisation models for the emission 
lines of Q1131+16. The ionisation parameter was varied over 

^ The version of CLOUDY used in this investigation is c08.00. 
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Figure 15. Diagnostic plot showing the measured [O 
III] (5007/4363) and [Fe vil] (6086/3759) ratios (filled point), com- 
pared with photoionisation model results for different densities 
and ionisation parameters. Equal densities are connected by the 
solid lines (as labelled), and the ionisation parameters are con- 
nected by the dotted lines (also labelled). 



the range -3.0 ^ log[U] ^ 0, in steps of log[U]=0.5, and the 
hydrogen density was varied over the range 3.0 ^ log(n£r 
cm""') ^ 8.0 in steps of log(n£f cm~^) =0.5. The rest of 
the properties of the model were kept constant: an ionising 
power-law index of -1.2 was assume^ll, and we only consid- 
ered radiation-bounded models that are optically thick to 
the ionisi ng continuum. In terms of abundances, Solar 84 wa s 
used (see lGrevesse fc Anderj|l989l . lGrevesse fc Noelslll993l '). 
Also no grains were included. This was so that the CLOUDY 
results were not affected by depletion of iron due to conden- 
sation onto dust grains. iNagao et all (|2003l ) have demon- 
strated for a large sample of AGN that the observed FHIL 
ratio^ cannot be reproduced in CLOUDY models that con- 
tain grains, but can be successfully reproduced with models 
not containing grains. 

The ratios measured for each of the main diagnostic ra- 
tios used in this study are given in Table 3. The WHT ratios 
are used in this investigation because the WHT spectrum 
has the greatest wavelength range. 

Considering first the low ionisation lines, we plot 
the transauroral [On](3726-f3729)/(7317+7330) and 
[Sn](6717-|-6731)/(4069+4076) ratios against each other 
in Figure 14. By comparison with the models, the ratios 
are consistent with hydrogen densities in the range 3.5 < 
log(n£f cm~^) < 4.5 at ionisation parameters -3.0 < log[U] 
< -2.0. The densities and ionisation parameters found 
for these low ionisation ratios are typical of these found 



^ Testing power-laws between -1.5 and -0.8, this power-law was 
the most successful at modelling the observed He IIA4686/H/3 
ratio compared to th e FHIL ratios . 

* The ratio used in iNagao et al.l l|2003l') was [Fe vii]A6086/[Ne 
V1A3426. 




Figure 16. Diagnostic plot showing the measured [Fe 
vil] (6086/3759) and [Fe vil] (5159/6086) ratios (filled point), com- 
pared with photoionisation model results for different densities 
and ionisation parameters. Equal densities are connected by the 
solid lines (as labelled), and the ionisation parameters are con- 
nected by the dotted lines (also labelled). 



Figure 17. Diagnostic plot showing the measured [O 
lll]A5007/H;3 and [Fe vil] (6086/3759) ratios (filled point), com- 
pared with photoionisation model results for different densities 
and ionisation parameters. Equal densities are connected by the 
solid lines (as labelled), and the ionisation parameters are con- 
nected by the dotted lines (also labelled). 



for the NLR of A GN in general (e.g. see iPetersonI Il997l . 
lOsterbrock fc Ferlan d 2006,). 

On the other hand, the high ionisation species 
show markedly different results. In Figure 15 the 
[Fo vii] (6086/3759) ratio is plotted against the [O 
III] (5007/4363) ratio. These ratios consistent with hydrogen 
densities 6.0 < log(nH cm~^) < 6.5 at ionisation param- 
eters in the range -1.5 < log[U] < 0. Similarly Figure 16 
shows the results found by comparing the [Fe vil] ratios of 
(6086/3759) and (5159/6086). These ratios consistent with 
hydrogen densities 5.5 < log(n// cm~"^) < 7.0 and ionisation 
parameters -1.0 < fog[U] < 0. The densities and ionisation 
parameters indicated by the latter ratios are much higher 
than those implied by the low ionisation emission lines. 

As mentioned previously, the [O in]A5007 emission line 
intensity is surprisingly low relative to H/3 for an AGN with 
a high ionisation emission line spectrum. Figure 17 illus- 
trates that such a low ratio is possible at high densities. The 
[O lll]A5007/H/3 emission ratio is consistent with hydrogen 
densities 6.5 < log(nfr cm~^) < 7.0 at ionisation parame- 
ters in the range -1.5 < log[U] < 0. At such densities, above 
the critical densities of the [O ni]AA5007,4959 transitions, 
the [O III] emission lines are suppressed relative to H/3; this 
explains the relatively low [O lll]A5007/H/3 ratio. 

FHILs with lower ionisation potentials, such as [Fe 
Vl]A5176, are also consistent with the results found using 
the [Fe vil] ratios. In Figure 18 the [Fe vil] (6086/3759) 
ratio is plotted against the [Fe vi]A5176/[Fe v]A3891 ratio. 
These ratios are consistent with hydrogen densities 6.5 < 
log(n£f cm~^) < 7.0 at ionisation parameters in the range 
-1.5 < log[U] < 0. 

Similarly, FHILs which are not from iron ions are con- 
sistent with the results found using the [Fe vil] ratios. Fig- 
ure 19 shows the [Ne v]A3426/[Fe vil]A6086 ratio plotted 



against the [Fe vil] (6086/3759) ratio. This diagram is con- 
sistent with hydrogen densities 5.5 < log(n£f cm~^) < 6.5 
and ionisation parameters -1.0 < log[U] < -0.5. Therefore 
the high densities and ionisation parameters found for the 
FHILs are not limited to iron ions. 

The HellA4686/H/3 ratio is relatively insensitive to 
the density of the emission region, but more sensitive 
to the shape of the ion ising continuum spectrum (e.g. 
IPenston fc Fosbu"rvl 1 19781 ). and whether the emi tting re- 
gions are matter bounded (jRobinson et al.l Il987h . There- 
fore, this ratio can be used to check whether the emission 
clouds are radiation bounded, and also whether the assumed 
ionising continuum shape is appropriate. The measured 
HellA4686/H/3 ratio is found to be 0.21±0.05 and 0.23±0.01 
for the WHT and Gemini spectra respectively. Figure 20 
shows this ratio plotted against the [Fe vil](6086/3759) ra- 
tio. The densities and ionisation parameters implied by this 
diagram are consistent with those deduced from the other di- 
agnostic diagrams involving the FHILs (5.0 < log(n_f/ cm~'^) 
< 7.0, -1.5 < log[U] < 0). In this case, the consistency be- 
tween the data and the models, as well as the results for 
other high ionisation line ratios, support the assumption of 
optically thick models with an ionising contirmum spectral 
index of -1.2. 

Although the models can successfully reproduce the 
FHIL line ratios in Table 3, the emission line ratios involv- 
ing lines with the highest critical densities (e.g. Ha/H/3, [Fe 
x]/[Fe vii]) cannot be successfully modelled using the range 
of physical conditions and ionisation parameters that suc- 
cessfully reproduce the [Fe vil] and [O ill] line ratios. Indeed 
the higher critical density line ratios can only be reproduced 
at higher densities {tih > cm~^) using CLOUDY. Figure 
21 shows Ha/H/3 ratio plotted against the [Fe x]A6374/[Fe 
vil]A6086 ratio. In this case, hydrogen densities of 7.5 < 
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Figure 18. Diagnostic plot showing the measured [Fe 
vil] (6086/3759) and [Fe vi]A5176/[Fe v]A3891 ratios (filled point), 
compared with photoionisation model results for different densi- 
ties and ionisation parameters. Equal densities are connected by 
the solid lines (as labelled), and the ionisation parameters are 
connected by the dotted lines (also labelled).. 

log(n/f cm~^) < 8.0 and ionisation parameters -1.5 < log[U] 
< -1.0 are required to reproduce the data. This suggests that 
some regions in Q1131+16 have even higher densities than 
implied by the [Fe vil], [Fe vi] and [O ill] diagnostic ratios. 

To summarise, most of the key high and low ionisation 
diagnostic ratios are consistent with the radiation-bounded 
photoionisation models. While the low ionisation lines re- 
quire low densities (3.5 < log(nH cm~^) < 4.5) and ionisa- 
tion parameters (-3.0 < log[U] < -2.0), the higher ionisation 
lines imply much higher densities (6.0 < log(nH cm"'') < 
8.0) and ionisation parameters (-1.5 < log[U] < 0). However, 
it is important to emphasise that our approach assumes that 
the lines contributing to a particular ratio are emitted by a 
single cloud with a single density and ionisation parame- 
ter; the reality is likely to be more complex. In particular, 
some of the low ionisation lines (e.g. [S ii]A4072, [Oii]A7330) 
may have contributions from the high density/high ionisa- 
tion regions, while some of the higher ionisation lines (e.g. 
[O lll]A5007, H/3) will have contributions from the low den- 
sity/low ionisation region. Therefore the differences deduced 
between the conditions in the high and low density regions 
may be even more extreme than implied by our simple anal- 
ysis. 



4.2 Explaining the high Ha/H/3 ratio 

The Ha/H/3 ratio gives an unusual result in Q113H-16. 
This may be due to the high density of the FHIL region 
in th is object coupled with a n extended partially ionised 
zone iGaskeU fc Ferlandlll984h . Such high densities can lead 
to significant coUisional excitation of the neutral hydro- 
gen atoms in the partially ionised zones of clouds in the 
FHIL region. This results in an enhanced Ha emission, but 
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Figure 19. Diagnostic plot showing the measured [Ne 
v]A3426/H;3 and [Fe vil] (6086/3759) ratios (filled point), com- 
pared with photoionisation model results for different densities 
and ionisation parameters. Equal densities are connected by the 
solid lines (as labelled), and the ionisation parameters are con- 
nected by the dotted lines (also labelled). 

no further enhancement in the other Balmer series emis- 
sion lines, which have lower cross-sectio ns for collision, as 
well as higher excitation temperatures (GaskeU fc Ferlandl 
iJi)84, Ostcrbrock & Fcrland 2006). Such a process could ex- 
plain the high Ha/H/? ratio measured in the spectrum of 
Q1131-I-16, since there is likely to be be a major contribu- 
tion from the FHIL region to the line flux of Ha. 

The Ha/H/3 flux ratio was investigated using CLOUDY, 
as described in ^ 4.1. The ratio indicated by the WHT data 
could only be reached at the highest densities {uh > 10^ 
cm"'': see Figure 21) and is not consistent with the con- 
ditions implied by the [Fe vil] and lower iron species ra- 
tios. Note that high Hq/H/? would also be favoured by lower 
metal abundances and a harder ionising continuum than we 
have assumed in our modelling (|Gaskell fc Ferlandl ITosi ). 

4.3 Location of the FHIL region 

The presence of FHILs in the spectra of the majority of ac- 
tive galaxies has led to a number of studies on their nature, 
kinem atics and location in AGN (e.g. see iMuUanev et al.l 
I2OO9I) . One of the most credible suggestions for the ori- 
gin of the FHIL emission i s th e inner torus wal l (see 
iMuravama fc Taniguchilll998l and iNagao et al.ll200ll ). The 
FHIL emission is rich in iron lines, which can be enhanced 
(relative to the classical NLR) by the evaporation of dust 
grains in t he inner torus wall , releasing the iron locked up in 
the grains (|Pier fc Voitll 19951 ). Moreover, studies of molecu- 
lar emission from the circum-nuclear molecular clouds of the 
centres of the Milky Way and other galaxies (hypothetically 
the torus) reveal densitie s in the range nH2 ~ 10^~® cm"'' 
(e.g. iPaglione et al.lll998l ). similiar to the densities deduced 
for the FHILs. 

In the case of Q1131+16, our results are consistent with 
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Figure 20. Diagnostic plot showing the measured He nA4686/H/9 
and [Fe vil] (6086/3759) ratios (filled point), compared with pho- 
toionisation model results for different densities and ionisation 
parameters. Equal densities are connected by the solid lines (as 
labelled), and the ionisation parameters are connected by the dot- 
ted lines (also labelled). 



Figure 21. Diagnostic plot showing the measured Ha/H/3 and 
[Fe x]A6374/[Fe vn]6086 ratios (filled point), compared with pho- 
toionisation model results for different densities and ionisation 
parameters. Equal densities are connected by the solid lines (as 
labelled), and the ionisation parameters are connected by the dot- 
ted lines (also labelled).. 



the idea that the FHILs are emitted by the far wall of the 
torus, observed with our hue of sight at a relatively large 
angle to the torus axis. We emphasise that such a geome- 
try is consistent with the relatively narrow line widths and 
small velocity shifts of the FHIL, since the circular velocities 
associated with the torus and any out-of-the-plane gas mo- 
tions would then be directed close to perpendicular to the 
line of sight. In this case, our results imply that the torus 
is dynamically cold, with velocity dispersion that is small 
relative to its circular velocity. 

It is difficult to entirely rule out the idea that the FHILs 
are emitted by an infalling molecular cloud, or a clump of 
the torus which has broken off from the rest of the torus 
and is falling towards the AGN. However, in such cases we 
might expect to see a larger velocity shift between the FHIL 
emission lines and the host galaxy rest-frame, unless the 
cloud happens to be falling perpendicular to the line of sight. 

Therefore the FHIL emission is most likely to origi- 
nate from the inner to rus wall of the AGN, supportin g the 
original suggestion of iMuravama fc Taniguchil l|l998l) and 



origmai suggestion 
iNagao et alT(|200ll ). 



4.4 The radial distance to the FHIL region 

The idea that the FHIL region is located in the torus can be 
further investigated by determining its radial distance from 
the AGN (rn). 

First we estimate vh based on the spatial information 
given in /3.1, and the parameters determined in ^4.1. The 
ionising luminosity (Lion) of the illuminating source is re- 
lated to the radial distance of the emission region from the 
ionising source (r) and the ionisation parameter (U), by: 



where nn is the hydrogen density, (/iz^);^^ is the mean ionis- 
ing photon energjij, and c is the speed of light. Then, solving 
for the radial distance of the low ionisation/density region, 
we get: 



r-L 



Lion 



47r [/i<c(M,„„ 



(2) 



and similarly for the high ionisation/density emission region: 



2 

TH = 



Li 



47r UHn»ic{hiy)^^^ 



(3) 



By combining equations 2 and 3, we obtain an equation 
which allows us to determine the relative radial distances of 
the high and low density regions from the AGN: 



(4) 



Lion = 'iirr UnHc{hu}ion; 



(1) 



Using the densities and ionisation parameters estimated in 
^4.1, we find that the radial distance of the FHIL emission 
region ranges from 0.0002 to 0.03 of that of the low ionisation 
region. As mentioned in ^3.6, the [O 11] A3727 emission is re- 
solved in the 2D spectrum of Q1131-I-16, and the half width 
at half maximum (HWHM) measured for the central peak of 
this emission line is found to be 0.53±0.05 arcseconds (cor- 
rected for the seeing). This corresponds to ri:, = 1.7±0.2 kpc, 
implying that the radial distance of the FHIL region from 
the AGN is in the range 0.30 < r_r/ < 50 pc — consistent 
with recent observational estimates of the size of the torus 
in nearby AGN, based on infrared interferometry and high 



^ We assume a value of 56 eV for {hy)^^^, which is based on an 
ionising power-law of -1.2, with photon energy limits of 13.6 eV 
and 5 keV. 
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resolution molecular line observations (e.g. IHaschick et al.l 
1 1994 IChou et al.ll2007l . iTristram et al.ll2007( ). However, we 
emphasise that this argument assumes that the [O ll] emis- 
sion line region is solely ionised by the AGN, which may 
not be true if stellar photoionisation is significant at a kpc 
scale. Also, if the contrast between the conditions in the low 
and high ionisation regions is larger than assumed here (see 
discussion in ^ 4.1), the radial distance of the FHIL region 
from the AGN will be correspondingly lower. Overall it is 
likely that this analysis, based on estimates of U and uh, 
overestimates the true radial distance of the FHIL region. 

If the ionising luminosity {Lion) is known it is also pos- 
sible to estimate the radial distance of the FHIL region from 
the AGN using equation 1 directly. Lion can be estimated 
from the bolometric luminosity {Lbol), which itself can 
be determined from the [O in]A5007 emission line luminos- 
ity (I / J OIII1 ) by manipulating equation 1 from iDicken et al.l 
1120091 1^ to give: 

T r 

Lbol = — i^[o\u\, \p) 

where C„ir is the NLR covering factor. Assuming typical 
NLR covering factor s in the range 0.02 < Cnlr < 0.08 
iNetzer fc Laoij|l993h . we estimate a bolometric luminosity 
in the range 2x10^* < Lbol < 10^* W. 

Alternatively, we can use our knowledge of the likely ge- 
ometry of the FHIL to estimate an upper limit on the cover- 
ing factor of the FHIL region and a lower limit on Lbol- As- 
suming that the FHILs are emitted by the inner wall of the 
torus on the far side of the AGN (see f 4.3), we see a max- 
imum of 25% of the total area of the inner wall of the torus 
(with the rest obscured). A typical torus with an opening 
angle of 45° has a total covering factor of Ctorus=0.7. There- 
fore, for our assumed geometry, an upper limit on the cov- 
ering factor of the FHIL region is Cfh/l =0.7x0.25=0. 175. 
Substituting this into equation 5, we obtain the following 
lower limit on the bolometric luminosity: Lbol > 10'^* W. 

We convert Lbo l to Lio n using the relationship given 
m lElvis et al.l [r994l {Lbol ~ 3.1 Lion), which leads to 
3x10^'^ < Lion < 3x10^* W. Then, using equation 1 and 
the model parameters in ^ 4.1, we estimate the radial dis- 
tance of the FHIL region from the AGN to be in the range 
0.3 < r£f < 30 pc. For comparison, based on our estimates 
of the AGN bolometric luminosity, the sublimation radius 
of the dust grains in the toruj^ is 0.26 < r^ub < 0.52 pc — 
consistent with the lower end of the range of our estimates 
of vh- This consistency supports the idea that the FHILs 
are emitted by the inner torus wall. 



4.5 The unusual strength of the FHIL 

It is also important to consider why the FHILs are so much 
stronger in Q1131+16 and other similar objects than in more 
typical AGN. 

First we must consider the possibility that, rather than 
an AGN, the spectrum of Q1131+16 is the result of a su- 
pernova, in particular a type Iln supernova. A type Iln su- 
pernova is a subclass of supernova for which it is believed 



that the supernova progenitor has undergone mass stripping 
of its outer layers in the late stages of its life. This pro- 
cess introduces a dense hydrogen gas in the circumstellar 
medium surrounding the progenitor. The Shockwave from 
the supernova interacts with the material, heating it. This 
results in a rich spectrum of narrow emission lines, which has 
some similar chara cteristics of the spectrum of Q1131-I-16 
(|Smith et al.l|20 0^). It is notable that early observations of 
the [ Fe vil] emission lin es in one such supernovae, SN 2005ip 
(see ISmith et al.1 120091') . indicate electron densities of 10 
cm"^ and higher (|Smith et al]|2009l '). 

One argument against this suggestion is that the 
[O ni]A5007 luminosity of Q1131-I-16 exceeds that of the 
hottest known typ e Iln supernova by 3 orders of magnitude 
(|Smith et al.1120091 ). Also the line emission from a supernova 
is expected to be vari able. In the c ase o f the type Iln su- 
pernova considered bv lSmith et~all (|2009l '). the [O in]A5007 
increased 8-fold after 100 days and then faded to half this 
intensity after 3 years. The [Fe vn]A5720 emission line de- 
creased 4-fold in intensity after 100 days, but after 3 years 
this intensity had increased to half its original intensity. As 
for the [Fe vn]A6086 emission line, this decreased in intensity 
by half after 100 days and then varied over 3 years around 
half its original intensity. However, in Q1131-I-16 there is no 
evidence for emission lines fading in the 3 years between the 
WHT and Gemini spectra. This lack of variability rules out 
any supernova component. 

The lack of variability in the emission lines of 
Q1131-I-16, also provides strong evidence against the the 
idea that its unusual FHILs are due to an illumination 
flare ca used by the t idal disruption of a star by the 
SMBH l|Komossa et al.1 12008 ). In the well-studied case of 
SDSSJ095209.56-H214313.3 jKomossa et al.1 l2009l i. FHILs 
such as [Fe vil]AA3759&5159, are seen to dramatically 
vary over a timescale of 3 years: the [Fe Vll]A3759 linqifl 
by 36% in intensity over 3 years, the [Fe vn]A5159 
line by 79%, and HenA4686 by 65%. In contrast, the 
same lines in Q1131+16 did not vary significant over 
a similar 3 year timescale. Moreover, the Balmer lines 
in Q1131-I-16 lack the variable, mult i-peaked profiles ob- 
se rved in SDSSJ095209 .56-K214313.3 (|Komossa et al.ll2008l 
fc lKomossa et al.ll2009l ). 

Rather than an illumination event, an explanation for 
the unusual strength of the FHILs may by provided by the 
inner torus wall being viewed at a specific orientation, at 
which the inner wall on the far side of the torus can be seen 
by the observer, but the quasar itself remains hidden, be- 
cause the torus is being viewed at a relatively large angle 
to its axis (see ^ 4.3). In this case, the visible area of the 
torus wall is high enough that the FHILs are strong, but 
the direct quasar continuum emission does not dilute the 
torus line emission (see Figure 22). Therefore the FHILs are 
observed with large equivalent widths. Viewed at smaller 
angles to the axis of the torus, the luminous quasar nuclear 
emission would become directly visible, substantially reduc- 
ing the equivalent widths of the FHILs. 

In this context we note that, despite the unusually large 



Assuming L[0jii] =5.51///^ rather than L[Qiii]=12L7f^. 
We calculate the sublimation radius using the expression for 
AGN sublimation radius given in lElitzudliool . 



Here all ratios are relative to [O lll]A5007. In 
SDSSJ095209.56-I-214313 .3 the [O ni] flux varies very little 
between the spectra (see lKomossa et al.ll2009l) . 
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Figure 22. (a) Schematic displaying an edge-on view of the possible observed orientation of the AGN of Q1131-I-16. (b) Schematic 
displaying the visible area of the inner torus wall of Q1131-I-16 as viewed by the observer. This special geometry will maximise the 
observed emission line flux of the inner wall of the torus yet still obscure the central quasar, thus lowering the underlying continuum, 
and hence resulting in FHILs with large equivalent widths. 



equivalent widths of its FHILs, the measured ratios of the 
FHILs to lower ionization lines (e.g. [Fe Vll]/[0 ill], [Fe 
vii]/[S ii], [Fe vii]/[0 i] and [Fe x]/[0 i]) for Q1131+16 faU 
well within the ranges measured for typical broad l ine AGN 

I Syl, Syl.5, NLSyl) in the sample compiled bv Nagao et al.l 
Therefore, viewed from a direction closer to the 
torus axis, we would expect Q1131-1-16 to display a nor- 
mal quasar/Seyfert 1 spectrum, but the multitude of fainter 
FHIL would be difficult to detect because of the dilution by 
the bright quasar continuum and broad line emission. 

However, it is more challenging to explain why the FHIL 
are significantly stronger relative to the low ioniza tion lines 
in Q 1131-H6 than in typical Type 2 AGN (Na gao et aH 
I2OO0I ). Given the geometry in Figure 22, we would expect 
to observe significant FHIL emission over a range of viewing 
angles, even if the quasar nucleus itself is obscured; rather 
than a sharp cut-off, for this simple geometry, the strength 
of the FHIL would gradually diminish as the orientation of 
the axis of the torus to the line of sight increased (i.e. the 
torus became more edge-on). Therefore, although the partic- 
ular geometry shown in Figure 22 would maximise the FHIL 
emission, it is unlikely that a specific geometry alone can ex- 
plain the rarity of objects like Q1131-I-16. We now discuss 
three further factors that might lead to the FHIL appearing 
relatively strong compared with the general population of 
Type 2 AGN. 

- A lack of conventional NLR emission. If the emis- 
sion from a "conventional" NLR on a scale of ~100 pc — 
1 kpc were strong, this would dilute the FHIL emission, and 
the overall spectrum would not appear unusual relative to 
other Type 2 AGN. Therefore, the conventional NLR emis- 
sion must be relatively weak in the case of Q1131-I-16 and 
other similar objects. Weak NLR emission can arise if the 
covering factor of the NLR is relatively small, or if the il- 
luminating AGN has only recently been triggered (within 
the last 100 - 1000 yr), so that the ionising radiation has 
not had time to traverse the scale of the NLR. However, 
while the lack of a conventional NLR might help to explain 
the unusual emission line ratios, it would not explain the 
unusually high equivalent widths measured for many of the 
emission lines in Q113H~16. Indeed, if this was the dominant 



factor we would expect the emission lines in Q1131-I-16 to 
have lower equivalent widths than most other Type 2 AGN, 
which is clearly not the case. 

- A lack of obscuration by large-scale dust struc- 
tures that are not coplanar with the torus. Along 
with a compact torus on a scale of 1 - 100 pc, many active 
galaxies show evidence for larger dust lane structures on a 
scale of 100 pc - 5 kpc. If the dust lanes are not exactly 
co-planar with the torus structures (e.g. thicker or different 
orientation), they will make it more difficult to observe the 
inner face of the torus, thus reducing the strength of the 
optical FHIL emission, and narrowing the range of angles 
over which the FHIL would appear strong. While the lack 
such large-scale dust structures would certainly increase the 
chance of observing unusually strong FHIL, their presence 
would in any case reduce the range of viewing angles over 
which the FHIL would be visible. 

- A recently triggered AGN. If the AGN has only re- 
cently been triggered then, not only might this result in re- 
duced emission from the conventional NLR (see above), but 
the inner wall of the torus itself might not yet have reached 
an equilibrium radius set by dust sublimation. In this case, 
the clouds close to the AGN are likely to be in t he process of 
evapo ration by the strong AGN radiation field (|Pier fc Voij 
Il995h and, depending on the density gradient, might have a 
higher ionisation parameter than clouds at the sublimation 
radius, leading to relatively strong FHIL. The problem with 
this explanation for Q1131-I-16 is that the radius of the FHIL 
emission line region is similar to the estimated sublimation 
radius of the torus (see f 4.4), suggesting an equilibrium 
situation. However, it is possible that future refinement of 
the physical conditions of the FHIL in Q1131-I-16 and the 
bolometric luminosity of its AGN might lead to a revision 
of this picture. 



Overall, while a specific viewing angle is likely to provide an 
important part of the explanation for the unusual spectrum 
of Q1131-I-16, one or more of the other factors discussed 
above may also be significant. 
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Table 4. T able comparing key properties determined for the FHIL regions in of Q 1131+16, III Zw 77 l|Osterbrocklll98ll 'l. Tololo 0109-383 
llFosburv fc Sansom 1983, .Muravaina k. Taniguchi.l99&') and ESO 138 Gl (.Alloin et al.lll992h . *shows systematic increase from low to 
high ionisation lines. 



Object 


Density 


Velocity Width 


[Oiii] (5007/4363) 




(cm~^) 


(FWHM, km s-l) 




Q1131+16 


10*>?iff>105 


361±29 


5.86±0.21 (WHT) 


III Zw 77 


ne>105 


~400* 


11.48±0.02 


Tololo 0109-383 


ne>10^ 


~330 


21.50±0.46 


ESO 138 Gl 


ne>105 


~1000 


30±2 



4.6 Q1131+16 in relation to other similar objects 

As mentioned previously, there are a few similar objects in 
which the FHILs are unusually strong (III Zw 77, Tololo 
0109-383 and ESO 138 Gl). The key properties of these ob- 
jects, as derived from their FHILs, are presented in Table 
4. In particular, we note that the densities derived from the 
FHILs are consistently high. Also all the objects show a par- 
ticularly low [O III](5007/4363) line ratio when compared to 
'typical' Seyfert galaxies (although in this respect Q1131+16 
is the most extreme). Finally, all the objects show relatively 
modest FHIL hne widths (FWHM), with the exception of 
ESO 138 Gl whose line widths are larger. 

Because there are relatively few objects with such un- 
usual features, it is entirely plausible that the properties 
of all these objects can be explained in a similar way e.g. 
with a specific viewing angle of the torus that maximises 
the observed emission area of the torus whilst still obscur- 
ingthe central engine, consistent with the scheme described 
by [Muravam a & Ta niguchi (^1998,) , .Nagao et al., (,2001, ') and 
iNagao et al] (j200ol ). 



5 CONCLUSIONS 

Q1131+16 has a remarkable spectrum which displays a mul- 
titude of FHILs of relatively high equivalent width, as well 
as the more common emission lines expected of an AGN 
spectrum. 

An in-depth study of the spectrum of Q113H-16 has 
revealed: 

• emission line ratios for the FHIL region that are con- 
sistent with high densities (10^ '' < uh < 10*" cm"^) and 
ionisation parameters (-1.5 < log[U] < 0); 

• similar velocity widths (FWHM) for the low and high 
ionisation lines, and a lack of evidence for substantial veloc- 
ity shifts of the FHILs relative to the galaxy rest frame; 

• a small radial distance from the AGN for the region 
emitting the FHILs (0.30 < rpHiL < 50 pc). 

Based on the high densities and relatively quiescent 
kinematics implied by our observations, it is likely that the 
FHILs in Q1131-I-16 are emitted by the torus wall, with the 
inner wall on the far side of the torus viewed directly by 
the observer, but the quasar itself remaining hidden. This 
geometry is also consistent with the relatively small radial 
distance found for the FHIL region. 

These results demonstrate the potential of the FHILs 
for probing the circum-nuclear obscuring regions in AGN. 
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Table 1: Line identifications for the nuclear spectra of Q1131+16. Flux 
ratios relative to H/3 are not corrected for reddening. Line IDs and red- 
shifts are determined via fitting a single Gaussian to the emission lines. 

Errors for the measured wavelengths and redsliifts include calibration er- 
rors as well as errors due to the measurement techniques. The total flux 
of the H/3 emission line is (4.38±0.11)xl0"^^ and (4.74±0.11)xl0"^^ ergs 
s^^ A^^ cm^'^ for the WHT and Gemini spectra respectively. Reference 
key for Table 5: (1)?, (2) ?, (3) Identified Using the FIVEL program, (4) 
Identified using NIST database, (5) ?, (6) ? & (7) ?. Notes key for Table 
1: (1) The broad line components were not modelled for these emission 
features to avoid negative fiux results for their line fiuxes. (2) The nar- 
row line components were not modelled for these emission features to 
avoid negative flux results for their line fiuxes. Notesl and Notes2 refer 
to the WHT and Gemini spectra respectively. 



Line ID 


A (A) 


Xwht (A) 


F/FHPwht 




zwht 




Xgemini (A) 


F/F H/3gemini 




zgemini 


Ref. 


Notesl 


Notes2 


[Mg V] 


2782.7 


3264.7 


± 


0.36 


0.158 


± 


0.024 


0. 


,1732 


± 


1, 


.lE-04 
























7 


- 


- 


Mg 11] 


2795.5 


3279.7 


± 


0.60 


0.630 


± 


0.126 


0. 


1732 


± 


1 


.5E-04 
























2 


- 


- 


Mg II] 


2802.7 


3287.9 


± 


0.71 


0.549 


± 


0.129 


0. 


1731 


± 


1, 


.5E-04 
























4 


- 


- 


[Fe IV] 


2829.3 


3319.6 


± 


0.60 


0.199 


± 


0.021 


0. 


1733 


± 


1 


.3E-04 
























4 


1 


- 


[Fe IV] 


2835.7 


3326.9 


± 


0.52 


0.419 


± 


0.021 


0. 


1732 


± 


1 


.OE-04 
























4 


1 


- 


O III 


3132.9 


3675.2 


± 


0.52 


0.358 


± 


0.027 


0. 


1731 


± 


1 


.4E-04 


3675.8 


± 


0.51 


0.382 


± 


0.041 


0. 


1733 


± 


1 


.4E-04 


1 


- 


- 


He I 


3187.7 


3740.1 


± 


0.61 


0.081 


± 


0.010 


0. 


1733 


± 


2 


.6E-04 


3738.7 


± 


0.70 


0.074 


± 


0.006 


0. 


1729 


± 


1 


.9E-04 


4 


2 


2 


He II 


3202.7 


3757.2 


± 


0.97 


0.095 


± 


0.024 


0. 


1731 


± 


2, 


.6E-04 


3757.6 


± 


0.55 


0.079 


± 


0.006 


0. 


1733 


± 


1, 


.5E-04 


6 


- 


- 


[0 HI] 


3312.3 


3886.6 


± 


0.64 


0.037 


± 


0.004 


0. 


1734 


± 


1 


.8E-04 


3885.0 


± 


2.10 


0.034 


± 


0.003 


0. 


1729 


± 


5 


.4E-04 


4 


1 


1 


[Ne V] 


3425.9 


4019.5 


± 


0.50 


1.010 


± 


0.085 


0. 


1733 


± 


1 


.2E-04 


4020.2 


± 


0.52 


1.051 


± 


0.035 


0. 


1735 


± 


1 


.3E-04 


1 






[N I] 


3466.0 
























4067.0 


± 


1.80 


0.026 


± 


0.003 


0. 


1734 


± 


4 


.4E-04 


6 




1 


He I 


3487.6 
























4092.3 


± 


1.70 


0.023 


± 


0.003 


0. 


1734 


± 


4 


.2E-04 


6 






[Fe VI] 


3662.5 


4297.2 


± 


0.51 


0.051 


± 


0.011 


0. 


1733 


± 


1, 


.4E-04 


4298.0 


± 


1.10 


0.042 


± 


0.006 


0. 


1735 


± 


2, 


.6E-04 


1 






HI 


3704.9 
























4346.7 


± 


0.52 


0.011 


± 


0.002 


0. 


1732 


± 


1 


.2E-04 


6 




2 


[O II] 


3726.0 


4371.9 


± 


0.61 


0.597 


± 


0.028 


0. 


1733 


± 


1 


.lE-04 


4371.6 


± 


0.59 


0.513 


± 


0.039 


0. 


1733 


± 


1 


.3E-04 


1 






[O II] 


3728.8 


4375.1 


± 


0.63 


0.156 


± 


0.049 


0. 


1733 


± 


1, 


.lE-04 


4374.7 


± 


0.60 


0.131 


± 


0.025 


0. 


,1732 


± 


1, 


.lE-04 


1 






[Fe VII] 


3758.9 


4410.5 


± 


0.51 


0.421 


± 


0.028 


0. 


1733 


± 


1 


.2E-04 


4410.9 


± 


0.51 


0.395 


± 


0.018 


0. 


1735 


± 


1, 


.2E-04 


1 






HI 


3797.9 
























4455.6 


± 


0.80 


0.036 


± 


0.004 


0. 


1732 


± 


1 


.8E-04 


6 






He I 


3819.3 
























4480.4 


± 


0.80 


0.023 


± 


0.003 


0. 


1731 


± 


1, 


.8E-04 


6 




2 


[Fe V] 


3839.3 


4504.1 


± 


0.91 


0.084 


± 


0.035 


0. 


1731 


± 


2 


.OE-04 


4502.6 


± 


0.52 


0.069 


± 


0.003 


0. 


1728 


± 


1 


.2E-04 


1 




2 


[Ne HI] 


3868.8 


4539.4 


± 


0.50 


1.103 


± 


0.050 


0. 


1733 


± 


1, 


.2E-04 


4539.3 


± 


0.79 


1.359 


± 


0.102 


0. 


1733 


± 


1, 


.7E-04 


1 






[Fe V] 


3891.3 


4565.6 


± 


0.56 


0.363 


± 


0.051 


0. 


1733 


± 


1 


.2E-04 


4565.4 


± 


0.60 


0.533 


± 


0.003 


0. 


1732 


± 


1 


.3E-04 


1 






[Ne HI] 


3967.5 


4655.3 


± 


0.51 


0.449 


± 


0.038 


0. 


1733 


± 


1 


.lE-04 


4655.7 


± 


0.59 


0.466 


± 


0.020 


0. 


,1735 


± 


1 


.3E-04 


1 






He I 


4026.2 
























4722.6 


± 


1.10 


0.013 


± 


0.002 


0. 


1730 


± 


2 


.3E-04 


6 




2 


[S II] 


4068.6 


4773.5 


± 


0.50 


0.098 


± 


0.012 


0. 


1733 


± 


1 


.2E-04 


4774.1 


± 


0.58 


0.142 


± 


0.012 


0. 


1734 


± 


1 


.2E-04 


1 






[S II] 


4076.3 


4782.5 


± 


0.50 


0.021 


± 


0.004 


0. 


1733 


± 


1, 


.OE-04 


4785.4 


± 


0.58 


0.019 


± 


0.006 


0. 


,1740 


± 


1, 


.2E-04 


3 


2 




m 


4101.7 


4812.7 


± 


0.51 


0.240 


± 


0.016 


0. 


1733 


± 


1 


.lE-04 


4812.3 


± 


0.50 


0.240 


± 


0.013 


0. 


1732 


± 


1 


.OE-04 


1 






[Fe V] 


4180.9 
























4904.1 


± 


0.52 


0.040 


± 


0.006 


0. 


1730 


± 


1 


.lE-04 


6 
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Table 1 - Continued 



to 



Line ID 


A (A) 


Xwht (A) 


F/F HiSwht 




zwht 




Xgemini (A) 


F/F H/3gemini 


zgemini 


Ref. 


Notesl 


Notes2 


[0 II] 


4317.2 










- 








- 






5066.6 ± 0.50 


0.056 ± 0.007 


0.1736 ± 9.9E-05 


6 


- 


- 


H7 


4340.5 


5092.6 


± 


0.50 


0.438 


± 


0.024 


0. 


1733 


± 


1 


.OE-04 


5092.7 ± 0.61 


0.454 ± 0.029 


0.1733 ± 1.2E-04 


1 


- 


- 


[0 III] 


4363.2 


5119.2 


± 


0.50 


0.899 


± 


0.019 


0. 


1733 


± 


1 


.OE-04 


5118.8 ± 0.59 


0.839 ± 0.031 


0.1732 ± 1.2E-04 


1 


- 


- 


[O II] 


4414.9 


5180.1 


± 


0.78 


0.067 


± 


0.023 


0. 


1733 


± 


1 


.OE-04 


5180.5 ± 0.55 


0.075 ± 0.011 


0.1734 ± l.lE-04 


4 


- 


- 


He I 


4471.4 










- 








- 






5246.4 ± 0.50 


0.044 ± 0.011 


0.1733 ± 9.5E-05 


6 


- 


- 


He II 


4685.7 


5497.5 


± 


0.63 


0.207 


± 


0.050 


0. 


1733 


± 


1 


.lE-04 


5497.9 ± 0.51 


0.234 ± 0.012 


0.1733 ± 9.3E-05 


1 


- 


- 


[Ar IV] 


4712.6 










_ 














5528.6 ± 0.67 


0.055 ± 0.010 


0.1732 ± 1.2E-04 


6 


- 


- 


H/3 


4861.3 


5703.5 


± 


0.50 


1.000 


± 


0.000 


0. 


1733 


■ 


1 


.OE-04 


5703.8 ± 0.52 


1.000 ± 0.000 


0.1733 ± 9.1E-05 


1 


- 


- 


[Fe VII] 


4893.4 


5740.8 


± 


0.51 


0.055 


± 


0.018 


0. 


1732 


± 


1 


.OE-04 


5740.6 ± 0.52 


0.057 ± 0.016 


0.1731 ± 1.2E-04 


4 


- 


- 


[Fe IV] 


4903.1 


5753.0 


± 


0.52 


0.012 


± 


0.004 


0. 


1733 


± 


1, 


.2E-04 


5753.1 ± 0.51 


0.038 ± 0.011 


0.1734 ± 1.2E-04 


4 


1 


1 


[0 HI] 


4958.9 


5817.7 


± 


0.50 


1.757 


± 


0.085 


0. 


1732 


± 


1 


.OE-04 


5818.1 ± 0.55 


1.914 ± 0.053 


0.1733 ± 9.5E-05 


1 


- 


- 


[O HI] 


5006.8 


5874.0 


± 


0.50 


5.270 


± 


0.155 


0. 


1732 


± 


1 


.OE-04 


5874.3 ± 0.55 


5.740 ± 0.158 


0.1733 ± 9.4E-05 


1 


- 


- 


[Fe VI] 


5145.8 


6036.9 


± 


0.52 


0.053 


± 


0.005 


0. 


1732 


± 


1 


.OE-04 


6036.7 ± 0.65 


0.047 ± 0.005 


0.1731 ± l.lE-04 


1 


- 


- 


[Fe VII] 


5159.0 


6052.7 


± 


0.81 


0.103 


± 


0.030 


0. 


1732 


± 


1 


.2E-04 


6052.7 ± 0.51 


0.123 ± 0.004 


0.1732 ± 8.4E-05 


1 


- 


- 


[Fe VI] 


5176.4 


6073.4 


± 


0.66 


0.097 


± 


0.024 


0. 


1733 


± 


1, 


.2E-04 


6072.9 ± 0.60 


0.094 ± 0.004 


0.1732 ± 9.9E-05 


1 


- 


- 


[Fe IV] 


5236.1 


6142.3 


± 


1.21 


0.035 


± 


0.005 


0. 


1731 


± 


2 


.lE-04 


6143.3 ± 0.62 


0.051 ± 0.006 


0.1733 ± 1. OE-04 


4 


1 


- 


[Fe VII] 


5276.4 


6189.6 


± 


0.65 


0.121 


± 


0.021 


0. 


1731 


± 


1 


.lE-04 


6189.3 ± 0.52 


0.113 ± 0.006 


0.1730 ± 8.4E-05 


4 


- 


- 


[Ca V] 


5309.1 


6228.6 


± 


0.70 


0.135 


± 


0.016 


0. 


1732 


± 


1, 


.lE-04 


6227.4 ± 2.10 


0.105 ± 0.009 


0.1730 ± 3.4E-04 


1 


- 


- 


[Fe VI] 


5335.2 


6259.9 


± 


2.25 


0.032 


± 


0.005 


0. 


1733 


± 


3 


.OE-04 


6259.7 ± 0.64 


0.046 ± 0.007 


0.1733 ± 1. OE-04 


4 


- 


- 


[Fe VI] 


5425.7 










- 








- 






6366.4 ± 0.64 


0.030 ± 0.004 


0.1734 ± 1. OE-04 


6 


- 


- 


[Fe VII] 


5720.7 


6712.0 


± 


0.56 


0.318 


± 


0.027 


0. 


1733 


± 


1 


.OE-04 


- 


- 


- 


1 


- 


- 


He I 


5875.6 


6893.1 


± 


0.66 


0.141 


± 


0.026 


0. 


1732 


± 


1 


.OE-04 


- 


- 


- 


1 


- 


- 


[Fe VII] 


6086.9 


7141.1 


± 


0.52 


0.495 


± 


0.028 


0. 


1732 


± 


1, 


.OE-04 


- 


- 


- 


1 


- 


- 


[0 I] 


6300.3 


7391.5 


± 


0.55 


0.196 


± 


0.017 


0. 


1732 


± 


1, 


.OE-04 


- 


- 


- 


1 


- 


- 


[S HI] 


6312.1 


7405.7 


± 


0.55 


0.060 


± 


0.016 


0. 


1733 


± 


1 


.OE-04 


- 


- 


- 


1 


- 


- 


[O I] 


6363.8 


7466.5 


± 


0.55 


0.063 


± 


0.006 


0. 


1733 


± 


1, 


.OE-04 


- 


- 


- 


1 


- 


- 


[FeX] 


6374.6 


7479.2 


± 


0.55 


0.234 


± 


0.034 


0. 


1733 


± 


1 


.OE-04 


- 


- 


- 


1 


- 


- 


[Nil] 


6548.1 


7682.2 


± 


0.61 


0.262 


± 


0.020 


0. 


1732 


± 


1, 


.8E-04 


- 


- 


- 


1 


- 


- 


Ha 


6562.8 


7699.5 


± 


0.70 


4.997 


± 


0.173 


0. 


1732 


± 


1 


.9E-04 


- 


- 


- 


1 


- 


- 




6562.8 


7663.4 ± 


16.06 


1.088 


± 


0.188 


0. 


1677 


± 


2 


.lE-03 


- 


- 


- 


1 


- 


- 


[N II] 


6583.4 


7723.6 


± 0.62 


0.786 


± 


0.059 


0. 


1732 


± 


1 


,8E-04 


- 


- 


- 


1 


- 


- 


[S II] 


6716.4 


7879.7 


± 0.65 


0.332 


± 


0.039 


0. 


1732 


± 


1 


.OE-04 


- 


- 


- 


1 


- 


- 


[S II] 


6730.8 


7896.1 


± 0.65 


0.219 


± 


0.039 


0. 


1732 


± 


1, 


.OE-04 


- 


- 


- 


1 


- 


- 


He I 


7065.7 


8290.0 


± 1.06 


0.111 


± 


0.024 


0. 


1733 


± 


1 


.3E-04 








4 






[Ar HI] 


7135.8 


8371.0 


± 0.69 


0.091 


± 


0.031 


0. 


1731 


± 


1 


.OE-04 








5 






[O II] 


7319.9 


8586.3 


± 0.99 


0.096 


± 


0.028 


0. 


1730 


± 


1 


.3E-04 








1 






[O II] 


7330.2 


8598.3 


± 0.99 


0.037 


± 


0.011 


0. 


1730 


± 


1 


.3E-04 








1 






[Ar I] 


7868.2 


9228.6 


± 1.03 


0.090 


± 


0.022 


0. 


1730 


± 


1, 


.3E-04 








4 






[Fe XI] 


7891.8 


9256.8 


± 0.67 


0.157 


± 


0.021 


0. 


1730 


± 


1 


.OE-04 








1 







o 
o 
o 

o 
o 



Table 1: Unidentified emission lines across both spectra. The flux ra- 
tios are relative to H/3. All presented wavelengths are measured in the 
observed frame. A notes key is included in Table 5. 

Xwht (A) Xgemini (A) F/F Hl3wht F/F HPgemini Notes Suggested ID 



3344.0±1.09 




0.141 ± 0.034 








2 


[Ar iv] 


3488.9 


± 





.91 




0.055 ± 0.008 








2 


[Ne v] 


3661.9 


± 





.52 


3662.9 ± 0.61 


0.072 ± 0.027 


0.091 


± 


0.007 


2 


[0 III] 










3747.9 ±1.10 




0.079 


± 


0.009 


. 


? 


3924.7 


± 


0, 


.51 


3925.1 ± 0.52 


0.420 ± 0.016 


0.434 


± 


0.034 




iii/[Ne v] 










3963.2 ± 0.56 




0.023 


± 


0.004 


2 


? 


4040.4 


± 





.66 


4039.4 ± 1.50 


0.122 ± 0.033 


0.099 


± 


0.036 


. 


O III 










4108.1 ± 1.50 




0.029 


± 


0.004 


1 


? 










4174.6 ± 0.60 




0.017 


± 


0.002 


2 


? 


4208.3 


± 


0, 


.51 


4208.9 ± 0.58 


0.164 ± 0.017 


0.242 


± 


0.012 




[Fe vii] 










4326.5 ± 0.52 




0.006 


± 


0.002 


2 


H I 










4336.2 ± 0.52 




0.008 


± 


0.002 


2 


H I 


4426.3 


± 


1, 


.24 


4425.3 ± 1.45 


0.010 ± 0.002 


0.008 


± 


0.001 




H I 










4438.4 ± 1.40 




0.017 


± 


0.004 


2 


? 










4863.6 ± 0.60 




0.034 


± 


0.006 




He I 










4904.1 ± 0.52 




0.040 


± 


0.006 




[Fe v] 










4933.8 ± 0.70 




0.030 


± 


0.006 




? 


4979.8 


± 


0, 


.57 


4979.7 ± 0.55 


0.091 ± 0.022 


0.080 


± 


0.009 




? 


5030.2 


± 


1, 


.13 


5029.6 ± 0.66 


0.072 ± 0.027 


0.047 


± 


0.014 




? 










5066.6 ± 0.50 




0.056 


± 


0.007 




[O II] 










5291.0 ± 2.40 




0.009 


± 


0.002 




[N in] 


5301.8 


± 





.59 


5302.1 ± 1.60 


0.056 ± 0.027 


0.030 


± 


0.004 


1 


? 










5443.5 ± 1.10 




0.024 


± 


0.006 




[N m] 










5544.5 ± 0.76 


- 


0.048 


± 


0.009 


- 


? 










5561.1 ± 0.63 




0.045 


± 


0.009 




? 










5649.8 ± 0.52 




0.034 


± 


0.006 




? 










5915.6 ± 0.60 




0.025 


± 


0.006 


1 


[Si II] 










6122.5 ± 0.60 




0.010 


± 


0.002 


2 


? 










6171.4 ± 0.81 




0.018 


± 


0.004 




? 










6238.3 ± 5.90 




0.025 


± 


0.007 




? 










6309.3 ± 0.78 




0.029 


± 


0.006 




? 










6350.1 ± 0.52 




0.036 


± 


0.004 




He II 


6481.8 


± 


1, 


.42 




0.023 ± 0.003 










? 


6495.8 


± 


1, 


.15 




0.044 ± 0.019 










? 


6588.3 


± 


2, 


.17 




0.024 ± 0.008 










? 


6605.3 


± 


2 


.18 




0.239 ± 0.016 










[Fe vi] 


6748.5 


± 


2, 


.94 




0.046 ± 0.010 








2 


[Nil] 


8403.8 


± 


3, 


.14 




0.062 ± 0.013 










? 



© 0000 RAS, MNRAS 000, 000-000 



